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The present work studies the quasi one-dimensional Ni(dmii) 2-hased 
compounds within a correlated model. More specifically, we focus our 
attention on the composed influence of the electronic dimerization- 
factor and the repulsion, on the transport properties and the local- 
ization of the electronic density in the ground-state. Those prop- 
erties are studied through the computation of the charge gaps (dif- 
ference between the ionization potential and the electro-affinity: IP- 
EA) and the long- and short-bond orders of an infinite quarter-filled 
chain within a t — J(t, U) model. The comparison between the com- 
puted gaps and the experimental activation energy of the semiconduc- 
tor NH2Me2[Ni(dmit)2]2 allows us to estimate the on-site electronic 
repulsion of the Ni(dmit)2 molecule to 1.16eV. 



I. INTRODUCTION 

Since the discovery of the Bechgaard salts in the 80 'si, the interest for such 
charge transfer organic solids is increasing steadily. The molecular character of 
these systems allows a very large diversity in their structural arrangement as well 
as in their chemical composition and leads to isery rich jphase diagrams exhibiting 
a wide variety of attractive physicaLpropertiesa (chargeEru and spinH,density wave 
phenomena, electronicta or structucalQ localizations, metal to insulatoicru or metal to 
superconductor phases transitionsa, etc . . . ). Among the latter, those related to the 
anisotropy of the conductivity (lXL_pseudo-lD or 2D) and hence, to the electronic 
instabilities of the Eermi surfaceBiij or to the strongly correlated character of the 
electronic structured receive considerable attention. The physical origin of the 
conduction band constitutes an important characteristic that classifies the different 
compounds into : 

• one-band 3 / 4-filled systems for which the conductivity is insured by the High- 
est Occupied Molecular Orbital (HOMO) of the donors — such as the TTF- 
based moleculesoo, 

• one-band 1 /4-filled systems for which the conductivity is insured by the Lowest 
Unoccupied Moleculap-Orbital (LUMO) of the acceptors — such as M(dmii)i 
and related moleculesE-3, 

• multi-bands systems for which the conductivity is insured both by the HOMO 
of the donor molecules and the LUMO of the acceptor oneslla. 
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These systems have been extensively and fruitfully studied using non correlated 
models sjich-as tight-binding (one-band model) and extended Hiickel (orbital model) 
methodsEO'Ea for periodic systems. The .topology of the inter-molecular interactions 
and the occurrence of a nesting vectoita on the Fermi surface is sufficient to un- 
derstand the. electronic instabilities that induce first order metal to insulator phase 
transitionsEErtj, such as Charge Density Waves and Peierls distortionso. However 
second order phase transitions, that see the metallic properties being gradually de- 
stroyed by electron-localization, cannot be explained without the introduction of 
electronic repulsionQ. More generally the interaction between the on-site electron 
repulsion U, the dimerization strength S and the electronic localization is a crucial 
problem for second order phase transitions that requires the treatment of the infinite 
system within a correlated model. The problem is further complicated by the diffi- 
culties in determining the numerical values of the repulsion integrals U. Indeed, these 
quantities are difficult to obtain from ab-initio calculations on fragments, since they 
should take into account, in an effective way, the effects of the long-range electron 
repulsions, whenever those are not explicitly treated in the model. This problem 
is well known since the 50 's by users of semi-empirical methods and an illustrative 
example is the drop of the effective Hubbard on-site tt- repulsion of carbon, atoms in 
conjugated molecules from lleV for the benzenetj to 6eV for polyenesE!J. An ex- 
perimental evaluation of the repulsion U is possible from the activation energies of 
insulating or semiconducting compounds. However this would require precise the- 
oretical evaluations of the charge gaps (ionization potential minus electro-affinity) 
as a function of U and the structural parameters. 

The physical properties of these charge transfer salts can be modelized within the 
Hubbard^ or t — J123 formalisms reduced to the conduction band. The resolution of 
even such simple Hamiltonians requires a method able to deal with infinite systems. 
The Density Matrix Renormalization Group (DMRG) methodEj provides such a 
procedure for one-dimensional problems. The aim of the present work will be to 
modelize acceptor-based one-dimensional charge transfer salts such as the isostrue- 
tural systems [Ni{dmit) 2 } 2 N R 4 where NR 4 = NMe 4 , NHMe 3 , NH 2 Me 2 . ■ xA. 
We will study a quarter-filled dimerized chain as a function of both the on-site re- 
pulsion U and the dimerization factor S. We will compute charge gaps and bond 
orders of this system, within an effective, {/-dependent, t — J model. Properties 
for donor-based 3 / 4-filled systems can easily be deduced by par±iple-hole symmetry. 
Since the NH 2 Me 2 [Ni{dmit) 2 ] 2 compound is a semiconductors we will be able to 
estimate the value of U for the Ni{dmii) 2 molecule involved in that type of infinite 
molecular crystal. 

In the next section we describe the DMRG method as well as the model Hamilto- 
nian. Section [II reports and analyses the results. The last section will be devoted 
to conclusions and perspectives. 



II. MODEL AND METHOD 



A. The model 



Single-band quarter-filled systems have an average occupation per site of only 
half an electron. Therefore, independently of the value of the on-site repulsion, it 
can be expected that the bivalent anion states associated to [Ni(dmit) 2 ] 2 ~ sites are 
very weakly represented in the ground-state wave-function. In the non-correlated 
limit (where the weight of the dianionic states is known to be maximum) their 
contribution is only 1/16 — 0.0625. It is therefore reasonable to eliminate explicit 
reference to the di-anions in the effective model and to allow only the neutral state 
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|0) (no electron on the LUMO), and the two singly anionic states | f) and | J,) (one 
electron on the LUMO) as accessible Valence-Bond states on each site. An often 
encountered confusion should however be clarified : the absence of the di-anions 
in the wave-function is not synonymous of an infinite on-site repulsion, as usually 
assumed by extension from the half-filled band case, but is only the consequence of 
their low statistical weight in a weakly-filled band system. In this case, the main 
physical effect of the bivalent anion molecular states is to lower the effective energy 
of the local singlets with respect to the local triplets. This is taken into account 
within at — J Hamiltonian, where the role of the J exchange integral provides this 
effect. 

H t -j = t Y^, ( a L a J> + a L ai > CT ) - 3 YL ( a !r a ii ~ a a a lt)KT a J.i - a a a j,r) C 1 ) 

<i,j> it <i,3> 
<i,j> a <i,j> 

where aJ CT (resp. aj iCT ) is the creation (resp. annihilation) operator of an electron 
of spin a on the molecule i, so that a\ jot ^ — a\ ^oj ^ (resp. a. t -[-dj-.j, — a^a^f ) is 
the creation (resp. annihilation) operator of a local singlet on the bond < i, j >. t 
refers to the hopping integral between adjacent (< i,j >) molecules i and j while 
J is the effective exchange between them. 



1. The effective exchange integral 

There are different possible ways to extract the value of the effective J as a 
function of the hopping and repulsion integrals. We choose to extract it in a way 
such as to reproduce the spectroscopy of a dimer in the Hubbard models. This 
choice is motivated by the fact that the model should be valid over the whole 
range of U/\t\, which excludes all perturbative techniques. The ground-state of the 
Hubbard Hamiltonian on a dimer is a symmetric singlet 1 E S : 

9 = — 71 — — 71 — ^ ' 



of energy EhubC^g 



u - Vu 2 + m 2 



while the first excited state is a triplet 3 E S of energy 0. Within the t — J model the 
first excited state is similarily a triplet of energy and the ground-state a symmetric 
singlet associated to the energy — 2 J. The equation of the low energy spectrum in 
the two formulations yields immediately 



Let us notice that the effective exchange integral is a function of the correlation ratio 
U/\t\, going from t in the uncorrelated limit to —2t 2 /U in the strongly correlated 
one. Even though J is extracted from the spectroscopy of a dimer, the validity of 
this model on infinite chains has been extensively studied in a previous workcJ. 
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2. The correlated effective bond order operator 



It is clear that the elimination of the explicit reference to the di-anions in the wave- 
function will obviously affect the values of any observable one may be interested in. 
One must therefore derive (as done above for the J) effective observable operators 
in order to take into account the effects of these states. 

In this paper we are more specifically interested in the bond order. It is easy 
to see that the bond order (over the < i,j > bond) has a diagonal matrix in the 
above S 9 , 3 S S basis set and that its value for the triplet is strictly zero, both in 
the Hubbard and t — J models. The value of the ground-state bond order in the 
Hubbard model is 

where 8 is defined in equation ||. This value will be assign to the singlet ground-state 
in the t — J representation. It comes an effective bond order operator 



P 2 

<i,j> o 



t 2 + j 2 

<hj> 



B. The DMRG method 

The DMRG is a very powerful method proposed by S. WhiteEl a few years ago for 
the treatment of one dimensional spin systems. Since then it has become one of the 
leading numerical tools for the study of quasi-lD correlated quantum systems. This 
success is due to both its excellent accuracy for systems as large as a few hundred of 
sites, and its flexibility in terms of the model (Heisenberg, t — J, Hubbard, Kondo, 
etc). 

In the DMRG approach the properties of the infinite system are derived by ex- 
trapolating the results of a succession of calculations on finite systems. Each one 
of these finite-system calculations is considered as a renormalization group (RG) 
iteration. The length N S i te of the chain increases very slowly at each iteration. For 
instance, in the present case, N s ite is increased by 2 sites (see Fig. |l|). In spite of 
the fact that N s u e increases, the dimension of the many-body Hilbert space is kept 
constant by means of the following procedure illustrated Fig. 0. 



m states 
n sites 



X 



X 



m states 
n sites 



m states 
n+1 sites 



X = X 



m states 
n+1 sites 



FIG. 1. Illustration of the renormalization procedure. A cross represents a Ni(dmit)2 
molecule or t—J site. Single (resp. double) lines represent weak (resp. strong) interactions. 
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At each RG iteration the system is divided in two side blocks spanned by m local 
states (typically, m = 60,82 and 100 in the present calculations) and two central 
sites which require mi local states, each to be represented exactly (e.g., mi = 4 
for a Hubbard model and mi = 3 for a t — J one). The Hilbert space for the 
N th iteration (i.e. for the iV s , te = 2N + 2 system) is obtained as the set of the 
antisymmetrized direct products of the four-blocks local states 

$ijki =\<Pi® $ ® <t>l ® 4>t®) 

where 4>\ is the i th local state of the block b. Thus the size of the Hilbert space is 
kept as (m x mi) 2 all along the renormalization scheme. 

Five main steps are involved in going from iteration N to iteration N + 1: 

1. the ground state of the N S i te = 2N + 2 system is obtained using a full config- 
urations interaction procedure, 

2. the reduced ground-state density matrix of the superblock formed by one of 
the side blocks and its neighboring site is calculated (see Fig. |l|), 

3. the reduced density matrix is diagonalized and the eigenvectors yielding the 
m largest eigenvalues (i.e., occupations) are obtained, 

4. the superblock local states space of size m x mi is projected onto the m most 
populated states derived in the previous step; the renormalized interactions 
within these superblocks, and between them and the central sites, are obtained 
by performing the corresponding unitary transformations, 

5. the resulting Hamiltonian is then used for calculating the ground-state at 
the following iteration, i.e., we return to the first step until convergence is 
achieved. 

The number of RG iterations performed in this work is 100, that is 202 effective 
Ni(dmit)2 molecules in the chain. The number of block states kept is 60, 82 and 100. 
The infinite system properties are then extrapolated from those of the finite (but 
large) systems. 

III. RESULTS AND DISCUSSION 

The M(drmt) 2-based systems do not present structural dimerizations. How- 
ever, molecular extended Huckel (EHT) calculations exhibit electronic dimeriza- 
tions through alternated hopping integrals along the chaincil. The recent reason- 
able description of the transport properties of the N H-iMe^\NiUimit)-^\'i and the 
N H Mez[Ni{dmit)2\2 systems by Fermi surface EHT calculations^ has proved the 
good quality of the calculated electronic dimerization ratio 5 = (t c — ti)/t c where 
t c refers to the mire-dimer hopping integral and t\ to the mier-dimer one. 8 — 
and 6=1 correspond respectively to the non-dimerized and to the product of 
dimers limits. As announced previously, we computed the charge gap A c , that is 
the extrapolated difference (toward the infinite system limit) between the ionization 
potential and the electron affinity as a function of both U/\t c \ and 5. If E(N S i te , N e ) 
is the energy of a finite N s u e sites system with N e electrons : 

A c = lim A c {N slte ) 

N 3 ite >oo 

= lim {E(N site ;N site /2 + 1) + E(N site ; N slte /2 - 1) - 2 * E(N slte ; N slte /2)} 
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. U/\t c \ spans the whole range of correlation strength between and +00, while S 
is varied between the two extreme values given by the real systems, that is S = 0.05 
for the metallic phase N H Me-i[Ni(dmit) 2)2 and 6 = 0.75 for the semiconducting 
one N H2Me2[Ni(dmit)2\2- In addition to the macroscopic picture given by the 
charge gap, we aimed to have a closer, microscopic understanding of the ground- 
state wave-function and of its localized versus delocalized character, that cannot 
be provided by extended Huckel calculations. These properties can be traced using 
the chemical bonding on the dimers, and more specifically the ratio between long- 
and short-bond order parameters 

, Pi v Pl(Nsue) 

A = — = lim — — r- 

Pc N M a—*co Pc{N S ite) 

where p c and pi are respectively the short- and long-bond orders. A is expected to 
vary from for the product of dimers limit to 1 for the totally delocalized system. 

Owing to the present renormalization scheme that increases the number of sites 
by only 2 at each iteration, one encounters two specificities. 

• The number of electrons is increased by 1 at each renormalization iteration, 
the ground state being a singlet for odd iterations (N = 2p — 1, Ap sites 
and 2p electrons) and a doublet for the even ones (N = 2p, Ap + 2 sites and 
2p+ 1) electrons. As a consequence all computed obscrvablcs present an even- 
odd iterations alternation, requiring to compute separately the odd and even 
iterations limits — even though they both should converge to the same value. 

• In order to keep a physically meaninfull system, one must impose the edge 
bonds to be consistently short. Thus, depending on the parity of the iteration, 
this leads to a long/short alternation of the central bond, on which the bond 
orders are computed. 

Three series of calculations have been performed for rn = 60, 82, and 100 block 
states kept. Let us notice that the m = 82, and 100 calculations treat exactly the 
finite systems up to 10 sites while the m = 60 only treats them exactly up to 8 sites. 

A. Charge Gap 

Figures ||a-b present the 5 dependence of the charge gap for 8 equidistant U/(U + 
At c ) values in the range [0, 1]. They correspond respectively to m = 60 and m = 100 
block states kept at each iteration. The results for m = 82 states have not been 
presented here since they are not qualitatively — and even quantitatively — much 
different from those obtained with m — 100. 
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60 states kept 




100 states kept 




FIG. 2. Charge gaps in t c units as a function of the dimerization factor 8. Each curve 
corresponds to a different value of J/t, therefore of p = U/(U + 4t c ) which varies between 
and 0.97 by step of 0.097. (a) m = 60 and (b) m = 100 block-states are kept. 

While the two figures exhibit quite similar results for large dimerizations and 
large U, they differ subtancially when either U or 6 are small. Indeed, in these 
latter regions, the m — 60 calculations largely overestimate the values of the gaps 
compared to the ones obtained with m — 100. In the large dimerization region 
- 5 > 0.5 — the increasing character of the gap as a function of the correlation 
strength U is well reproduced, eventhough the discrepancy with the m = 100 curves 
increases rapidly when U decreases. In the small dimerization limit however, the 
m = 60 results do not even exhibit the correct behavior as a function of U . In- 
deed, for <5 = 0.05, the smallest computed gap corresponds to the largest value of 
U (the extrapolated gap for U — and 5 = is 0.05 instead of 0). Moreover, 
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unphysical crossings between ?7-curves can be observed. Those artefactual effects 
can be directly related to an insufficient number of block states kept in the m = 60 
calculations, since they disappear when m is increased to 82 or 100. This larger m 
value required for an adequate treatment in the small S and/or small U region can 
be rationalized in the following way. 

• The populations of the local block states in the total wave-function arc 
knownc3 to follow a decreasing exponential of the form 7 = expL-ra(U)n), 
where a(U) is a decreasing function of U. As documented in refe3, the ac- 
curacy of a DMRG calculation is directly related to the fraction of the total 
population provided by the m block-states kept. Hence, as U decreases, a 
larger number of states m needs to be kept. 

• An optimal use of the DMRG procedure supposes that the perturbation im- 
posed on the electronic wave-function follows a linear dependence with the 
system size. In the present case, the addition or removal of an electron is a 
major perturbation for small size systems (50% of the electronic population 
in the first iteration for the N S i te /2 — 1 system) whereas it is a negligible one 
for larger systems (100/ (N + 1)% at the iV th iteration). This problem can 
only be overcome by insuring that all local states needed for the description 
of the large systems (weakly perturbed) are not removed in the first few itera- 
tions. That is by increasing the size of the largest system treated exactly and 
equivalently enlarging m (8 sites for m = 60, 10 sites for m — 82 or m — 100). 

From now we will focus our attention on the m = 100 results. 

As expected the charge gap is an increasing function of both the dimerization ratio 
and the on-site repulsion. In agreement with the quarter-filled dimerized Hiickel 
model, the U = curve does not exhibit any significant gap. Its maximum computed 
value is 0.005 |t c | and can be considered as characteristic of the calculations precision. 
Similarly, all curves extrapolate towards a zero gap for S — 0. 

In large U limit the exchange J parameter tends to ; the infinite-!/ Hamiltonian 
is therefore limited to the kinetic part acting on a Hilbert space excluding all doubly 
occupied sites. 

Hu^ + oo = ( a i<T a 3° + a )* a ™) 

<i,j> a 

The spin part of the latter Hamiltonian is completely uncoupled from the -space 
part and the system is equivalent to a spinless fermion one. As shown in ref.E3 the 
ground-state energy per site can be expressed as 

where 

1 ) 7rn for n > 1/2 

~ \ vr(l -n) for n < 1/2 

n refers to the band filling, E to the elliptic integrals of the second kind. One sees 
immediatly that the energy derivative presents a discontinuity as a function of the 
band-filling, discontinuity that can be associated to the charge gap. 
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Turning back our attention to fig. one sees that our calculations converge well 
to the above expected limit for U — > +00. 

One of the main interest of having an accurate estimate of the charge gap within 
a correlated model is the ability it provides to evaluate the effective repulsion U . 
Organic conductors are well-known to behave as strongly correlated systems, since 
the rather weak inter-molecular interactions lead to quite large U /\t\ ratio. However 
the difficulties encountered to obtain reliable U has led to very large controversies 
for yearsu. On one hand, a non biased extraction from experimental datas would 
require the use of a correlated interpretative model beyond the standard non or 
weakly correlated methods (Hiickel, RPAcj). On another hand, direct computa- 
tion from ab initio quantum chemistry methods on the molecular unit is expected 
to (largely) overestimate the value of U . As already mentionned in the introduc- 
tion, the U integral is an effective integral that takes into account number of other 
correlation effects than the purely one-site repulsion. 

The availability of accurate gap values within a correlated model, valid over the 
whole U/\t c \ range, allows us to confront them to the experimental activation ener- 
gies in order to extract reliable values of U. Single-crystal temperature-dependent 
conductivity measurements have disclosed an activation energy of E a — 0.21eVcj, 
at room temperature, for the semiconductor NHiMei \Ni(dmi£)<i\i- Considering 
that at room temperature the gap is reduced by the thermal energy, the fc&T value 
should be added to E a in order to obtain its T — estimation. One therefore 
obtains an estimated experimental gap at T = OK of E a (0K) = 0.22eV, that-is 
E a (0K) = 0.76|i c | (where t c = 0.29eU has been taken from EHT calculations^) . 
This corresponds to U = 4.0|i c | = 1.16eV for 5 = 0.75. Following the same line of 
argument for the NHMe^ [Ni{dmii)2\2 compound, the \k\jT factor appears to be 
larger than the computed gap for the whole range of U/\t c \. Although the calculated 
gap A c never equals zero, the small values it takes are in agreement with the metal- 
lic behavior exhibited by this system. Using now the value of U obtained for the 
Ni(dmit) 2 molecule and the EHT t c = 0.237eV value for the NHMe 3 [Ni(dmit) 2 }2 
system, the correlation ratio comes to be U = 4.89|t c |, slightly larger than for the 
semiconducting phase. 



B. Bond Order 



Fig. U reports the short- and long-bond orders (respectively referenced as p c and 
pi) as a function of U/\t c \ and for different dimerizations. p c and pi have been 



computed according to the expression derived in section II A 2, i.e. the effects of 
the di-anionic configurations are taken into account in an effective way. From a 
technical point of view, the bond order was computed on the central bond at each 
iteration. In the present rcnormalization scheme, this bond alternates from long for 
odd iterations, to short for even iterations. Additional oscillations occur between 
(4p)- and (4p+2)-electrons systems for long bonds and (4p+l)- and (4p+3)-electrons 
systems for short bonds. We therefore performed 4 independent extrapolations. For 
short bonds, the (4p + 1)- and {Ap + 3)-electrons extrapolations converge towards 
the same values all over the U/\t c \ range. However, as can be seen on Fig. ||, for 
long bonds, the (4p + 2)- and (4p)-electrons extrapolated_23/ curves split over for 
small U, following the aromatic, anti-aromatic alternationEl 

As expected, when the dimerization factor goes to 1, the system tends toward 
a product of independent singly-occupied dimers, that is p c tends to 0.5 and pi to 
0. Looking now at the U = and U — > +00 limits, one can perform an analytical 
calculation as a function of 5. It is easy to reach the Hiickel result 
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t io V 1 + (1 - <*) 2 + 2(1 - S) cos6 

The_J7 — ► +00 result can be obtained from the spinless fermion case as treated in 
ref.E3, and 



p c (S) = < eld > 
1 

~ 2tt 



±{5F(l/2,q) + (2-5)E(l/2,q)} (6) 



where 



«r/a 1 
E{12,q)= - . d9 

Jo Vl-4gsin 2 6» 

r /2 

E(l/2,q)= Jl - 4q sin 2 6 de 

Jo 
1-6 

q 



(2-sy 

The asymptotic values are represented by open circles in Fig. pi all of them in very 
good agreement with our computed curves. 

According to the dimerization factor S, the different systems studied in fig. || 
present two types of behavior. The first type, which will be referred as localized, 
includes the 8 = 0.5 and 8 = 0.75 (NH2Me2 [Ni(dmit) 2)2) systems. It is character- 
ized not only by a small pi/p c ratio (as seen in Fig. ^) but also by a non monotonic 
behavior of the short-bond order as a function of U/\t c \. Starting from the Hiickel 
limit, p c begins to decrease, goes through a minimum (p™ m (i5)), then increases up to 
a U — > +00 limit, slightly smaller than the U = one. The U/\t c \ value associated 
to p™ m (S) is a decreasing function of 8. The second type of behavior, to which we 
will refer as delocalized, is represented by the 8 = 0.05 (NHMez [Ni{dmit)2[2) sys- 
tem. It is characterized by a pi/p c ratio close to 1 and a monotonically decreasing 
p c curve as U/\t c \ increases. The 8 — 0.2 p c curve however exhibit a intermedi- 
ate behavior between the two previous ones, p c first decreases as U increases, then 
reaches an asymtptotic value around U = 4|£ c |. 
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FIG. 3. Long- and short-bond orders as a function of U/(U +4t c ). Solid line: short-bond 
order, dashed line: long-bond order. The circles correspond to analytical evaluations of 
the asymptotic values for (7 = and U — ► oo. 



The long-bond order present a qualitatively (but not quantitatively) similar be- 
havior for all dimerization values, decreasing as the correlation strength raises and 
exhibiting an inflexion point. The latter occurs for smaller U/\t c \ values as S in- 
creases. As expected the long-bond order diminishes as 5 increases, while the reverse 
tendancy is observed for the short-bond order. 

Figure || presents the pi/p c ratio as a function of the correlation strength. It 
can be seen in this figure that — as expected — the correlation has a localizing 
effect since pi/p c diminishes as £7/|f c | increases for all values of 5. It is more in- 
teresting to compare the hopping-integral ratio (horizontal lines) to the bond-order 
one. The delocalized system exhibit similar ti/t c and pi/p c ratios. On the contrary, 
the localized systems show much smaller long- to short-bond order ratios than the 
hopping-integrals do, i.e the collective effects seem to strongly enhance the local- 
ization of the electron density on the dimers. 



11 



* * NHMe 3 

^-^5=0.2 
* — * §=0.5 
* — *NH,Me, 



V '"0.0 0.2 0.4 0.6 0.8 1.0 
U/(U+4tJ 

FIG. 4. Bond-orders ratio: pi/p c - Solid line: 5 = 0.05, dash-dotted line: S = 0.20, 
dotted line: S = 0.50, long-dashed line: S — 0.75. 

IV. CONCLUSION 

The present work deals with the correlation and dimerization effects in quarter- 
filled quasi one-dimensional organic conductors such as the X [M(dmit)2\2 radical 
anions compounds (where M stands for M, Pt or Pd). We treated them as a 
correlated dimerized chain within a t — J model. The effective exchange integral 
J is taken as a function of the hopping integral t and the correlation strength U, 
derived to reproduce the dimer low energy spectroscopy. We computed accurately 
the charge gaps for the infinite chain (using the DMRG method) over the whole 
range of electronic repulsion and dimerization factor. The comparison between 
our theoretical results and the experimental activation energy on semiconducting 
Ni(dmit) 2 -based compounds allowed us to extract a reasonable evaluation of the 
one-site repulsion integral U for the Ni(dmit)2 molecule: 1.16ey, i.e. U/\t c \ = 4.0 
for the semiconducting NH^Me^ [Ni(dmii)2\2 compound, and U/\t c \ = 4.89 for 
the metallic NHMe^ [M(c?TO?t) 2 ]2 one. It would be interesting to extract, in the 
same way, the U/\t\ value for other molecules widely encountered in such molecular 
compounds, specifically since there is a large number of controversies about their 
numerical values. 

We also computed the long and short bond order as a function of both S and U 
and showed that the collective effects enhance, in strongly dimerized systems, the 
electronic localization on the dimer units. Meanwhile an unexpected C//|t c | behavior 
has be found for these dimerized systems, with the observation of a minimum on 
the short-bond order curves. The increase of p c for large U support the idea that 
the electron repulsion enhance the accumulation of electrons on the strong dimcrs, 
leading to a product of dimer zeroth-order picture even for moderately dimerized 
systems (<5 = 0.5). 
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